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Abstract 

Let K be an imaginary quadratic number field with class number 1. We 
describe a new, essentially linear-time algorithm, to list all isomorphism 
classes of cubic extensions L/K \xp to a. bound X on the norm of the 
relative discriminant ideal. The main tools are Taniguchi's |18) general- 
ization of Davenport-Heilbronn parametrisation of cubic extensions, and 
reduction theory for binary cubic forms over imaginary quadratic fields. 
Finally, we give numerical data for K = Q(i), and we compare our results 
with ray class field algorithm ones, and with asymptotic heuristics, based 
on a generalization of Roberts' conjecture [19] . 

1 Introduction 

Given a number field a positive integer n and X > 0, we define J^K,n{X) to 
be the set of isomorphism classes of extensions L/K such that 

[L:K]^ n, and Mk/q{KL/K)) < X, 

where X){L/K) is the relative discriminant ideal of the extension L/K. Sets 
of this type may be enumerated algorithmically (usually over Q) using the ge- 
ometry of numbers, following Hunter-Martinet's theorem |14| . Asymptotically, 
their cardinality as X tends to infinity is the subject of folklore conjectures, 
predicting for instance that it should be of the order of X , strikingly refined by 
Malle [121 who also fixes the Galois group of the Galois closure of L/K. Small 
values of n are of particular interest, since computer tests become comparatively 
easier and more theoretical results are available; see [2] for a recent survey. 

In the present paper, we will focus on the case n = 3. Belabas's algorithm [T] 
lists all representatives of J-q^3(X), in time Oe(X^+'^), essentially linear in the 
size of the output. We consider the problem of generalizing this algorithm to 
other base fields and we will solve it completely when K is imaginary quadratic, 
with class number 1. Our main result is as follows: 

Theorem. Let K he an imaginary quadratic number field with class number 
hji = 1. There exists an algorithm which lists all cubic extensions in J^K,3iX) 
in time Oe(X^+^), for all e > 0. 

For an arbitrary fixed number field K, Datskovsky and Wright [SI Theorem I.l] 
proved that the cardinality of J-k,3{X) is asymptotic to a constant (depending 
on K) times X as X — cx). It follows: 
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Corollary. The algorithm runs in time essentially linear in the size of the 
output. 

The algorithm uses two main ingredients : 1) a general description of isomor- 
phism classes of cubic extensions L/K as classes of suitable binary quadratic 
forms in K[x,y] modulo a GL2 action; 2) classical reduction theory in the spe- 
cial case where K is imaginary quadratic. Enumerating cubic extensions then 
amounts to enumerating integer points in an explicit fundamental domain, cut 
out by the extra condition MK/q{^{L/ K)) < X. 

It is interesting to compare our algorithm with the classical one, using class 
field theory (see Section 9.2.3 of [4J): the latter works in time O^iX^^^'^''), 
unless we assume the Generalized Riemann Hypothesis to obtain Oe(X^+^). So 
our algorithm has better unconditional complexity. Moreover, even assuming 
GRH, as we did in our PARI/GP implementation, the ray class field algorithm 
is slower than ours (see section [5]). 

Section[5]is devoted to our two ingredients : Taniguchi's theorem [T5], which 
generalizes the Davenport- Heilbronn bijection used by Belabas [IJ, and general 
facts about reduction theory for integral binary cubic forms over imaginary 
quadratic fields. In Section [3] we further assume that K has class number 1 
and study the action of GL2{Ok) on binary cubic forms and obtain a specific 
fundamental domain, as well as explicit numerical bounds for the coefficients 
of reduced forms. Section 0] describes the core of our algorithm and Section [5] 
explores in detail the technical issues encountered during the implementation 
of the algorithm. The final section [6] presents some timings for our PARI/GP 
implementation, over K = Q{i). 
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2 Notations and preliminary results 

In this section, we recall known results, needed for our algorithm. 
2.1 Taniguchi's theorem 

Definition 2.1. Let O be a Dedekind domain, and let K be its quotient field. 
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• Let C{0) be the set of "cubic algebras" that is, isomorphism classes of 
O-algebras that are projective of rank 3 as O-modules. 

• For every fractional ideal a of O we define 

C{0, a) = {Re C{0) I St(i?) = the ideal class of a}, 

where St(i?) G C'l{0) is the Steinitz class of R, thus R is of the form 
ujiO © UJ2O © waa, for appropriate wi, 012, 013 e Frac(i?) := R (g)© K. We 
define the discriminant ideal d{R) — disc(ti;i,a;2, ^3)0^, where as usual 
disc(wi, W2,a;3) = det TrFrac(i?.)/_R'(wiWj). 




7 e a 6 e O 



Va ^ {F {a,b,c,d) \ a e a,b e 0,c Q a \dea ^}. 

If F e Va, its discriminant disc(i^) = b"^ - 27 a"^ d"^ + ISabcd - iac^ ~ Ab^ d 
belongs to a"^ . 

• We consider elements of Va as binary cubic forms, under the identification 
(a, 6, c, d) — ax^ + bx^y + cxy^ + dy'^ and we define a left-action of Ga on 
K by 

M -F ^ (det My^F{ax + /3y, jx + Sy), 

where M = ^) *^ 

The following theorem generalizes the Davenport-Heilbronn [9J theory, corre- 
sponding to the special case O = Z, to cubic algebras over an arbitrary Dedekind 
domain O: 

Theorem 2.2 (Taniguchi [18]). There exists a canonical bijection between C{0, a) 
and Va/Ga such that the following diagram is commutative: 

K/G„ > C{0,a) 

discj^ ^ 

0-2/(0^)2 > {integral ideals of O} 

where d is the relative discriminant ideal map. 



Remarks. 

• A computation proves that the vertical "disc" is well defined. The other 
vertical map 5 is well-defined since an O-algebra isomorphism preserves 
the discriminant. 

• We slightly changed the notation from Taniguchi's paper, to keep consis- 
tent with the notation of the following sections (Taniguchi's action M * F 
is given by (M*) • F). 
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Corollary 2.3. Let K be a number field with class number hx — 1- Let O = 
Ok be its ring of integers. Then Taniguchi 's bijection simplifies to a bijection 
between binary cubic forms with coefficients in O modulo GIj2{0) and cubic 
O-algebras. 

To enumerate relative cubic extensions L/K, we shall select only the cubic 
O-algebras R which are both domains and integrally closed: those algebras are 
exactly the classes of the Ol- The algebra i? is a domain if and only if F is 
irreducible over K. Being integrally closed is a local property; it is equivalent 
to p-maximality at all prime ideals p C Ok such that | c)(i?) and this can be 
tested using Dedekind's criterion [4, Theorem 2.4.8]. As was done in [1], it is 
possible to use sieve methods to control the complexity of this step by avoiding 
costly discriminant factorizations. 



2.2 Fundamental domains in hyperbolic 3-space 



In this section, we describe fundamental domains for the action of Bianchi groups 
on hyperbolic 3-space, which underlie the reduction of binary Hermitian and 
cubic forms (to be dealt with in the next two sections). 

Definition 2.4. Let H = K + Mi + Rj + Rfc be the algebra of quaternions, let 
C = M + Ri be the subfield of complex numbers, and let 

Hz = {z + tj\z<E<C,t<EM:\} 

= {h = z + tj \ h Cz M, such that the k-component is 0,t > 0}, 

denote hyperbolic 3-space. We define the action o/PGL2(C) on by M ■ (z + 
tj) {z' + t'j), with 



where M 



A 

C D, 



p^AC + zAD + zBC + BD 

p^\C\^ + zCD + zCD + \D\^ 
\det{M)\t 

p'^\C\^ + zCD + zCD + \D\2' 
e PGL2(C) andp'^ = \z\^ + t^ . 



(1) 



t' = 



Remark. With the quaternion notations (and operations), this translates to 
the neater formula 

M ■h={Ah + B){Ch + D)-^. 

Definition 2.5. Let K — Q{^/dK) be an imaginary quadratic field of discrimi- 
nant dK < and class number 1. We define 



Fn 



'i) 



zeC 



Fn 



-3) 



zeC 



< Re{z) < ^,0 < Im(z) < ^ 



-l/2<Rc(z)<i 0<lm(z)<^|, 
< Re(z) < i Re(z) < Im(z) < ^ Re(z) 



Fk = |z e C |Re(z)| < 1/2,0 < Im(z) < , w/ien 7^ -2,-3,-4. 
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Moreover we set 



Bk = {z + tj eH3 



z G Fk and \z\ 



t^ > 1 



Let 



Fa. 



F n 



F, 



{z + t] e Bk I Re(z) < Im(z) i/ z + tj £ OBk } 
{z + tjeBxl Re{z)>Oifz + tjedBK} 
{z + tjeBxl lm{z) > if z + tj e dBx} , 



where BBk denotes the boundary of Bk- 

Finally for K such that dK ^ —3, —4, we define 



Fk 



z + tj e Bk 



Re(z) < P/4 if [lm{z) = and |Re(z)| < |z| 

Re(z) >0 tf (\z\^ + <2 = 1 or |Re(z)| = 1/2 or Im(z) = 



Theorem 2.6. Let K he an imaginary quadratic number field of class number 
1, let O he its maximal order, and let Fk be as defined above. 

1. Fk is a fundamental domain for the action of VG\j2{0) on Us- No two 
points in Fk are PGh2{0)-equivalent. 

2. There exists a constant tK > such that t > tK for every z + tj G Fk ■ 
The value of t\ is given in the following tables : 



3/4 



D 


1 


2 


3 


7 


11 




1/2 


1/4 


2/3 


3/7 


2/11 



D 19 


43 


67 


163 


2/19 


2/43 


2/67 


2/163 



Proof. 1. Since PGL2(e')/PSL2(0)~C/(C)2 our hypotheses imply that 
its cardinahty is 2. Using the weh- known fundamental domains for the 
PSL2(C') action on T-L-^ (sec for example [lOj) the result follows. 

The action of PGL2(0) on the boundary of is generated by the fol- 
lowing matrices: 



(a) Either 



-1 

1 



acting on points z + jt such that \z\ +t^ ~ 1. 



(b) Or translations of the form 



1 a 
1 



, for an appropriate a G C 



A tedious computation yields the result. 
2. See [S] and [20] for details. 



□ 
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Remark. Thanks to Definition 12.51 and Theorem 12.61 we have explicit bounds 
for z and i-components of elements in a fundamental domain of "Hs modulo 
GL2(C'), when O is principal. Unfortunately, when Kk ^ 1, we do not have a 
lower bound for t (there are points in the boundary of the fundamental domain 
such that t = 0), and this will prevent us from bounding the coefficients of 
reduced forms. This is the reason why we will restrict our work to the class 
number 1 case. 



2.3 Reduction of binary Hermitian forms 

Before tackling cubic forms, we recall the classical reduction theory of binary 
Hermitian forms modulo GL2(C'), where O is the maximal order of an imaginary 
quadratic field. 

Definition 2.7. Let {P,Q,R) denote the binary Hermitian form 

H{x, y) = P\x\^ + Qxy + Qxy + R\y\^ , P,Q,Re C, 

of discriminant disc(i/) := —A = jQp — PR, and let he the set of positive 
definite binary quadratic Hermitian forms over C; in other words, 

^ = {{P,Q,R): P,ReR+,Q e C, disc(P, Q, R) < O} . 

The group PGL2(C') acts on ^ via 

M ■ H{x,y) = H{ax + I3y,^x + 5y), where M = ^) e PGLzCO). 
Remark. It is customary to identify the Hermitian form 

<^ " (5 i) (:) 

with the Hermitian matrix H = (^^'^ ; the PGL2(C') action is then 

M ■H = M* X H X M, where M* = (M)*. 
Lemma 2.8. Let^ : ^ /R*^ Hj be defined by: 



$((F,Q,i?)) = -| + ^j. (2) 

^ is a bisection wich commutes with the action o/PGL2(C'). 

This defines natural representatives for orbits of Hermitian forms modulo 
PGL2(0). Namely 

Definition 2.9. Let H ^ a binary hermitian form. H is called reduced if 
and only if <^{H) G Tk. 

Lemma 2.10. Let {P, Q,R) — P \x\^ + Qxy + Qxy + R \y\^ be a reduced Her- 
mitian form in ^ , with discriminant — A = |Q| — PR- We have 



P<^. (3) 
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\Q? < ckP\ (4) 

and 

Pi?<(l + j)A, (5) 
where ck is a constant depending only on the number field K , defined as follows 



Ck 



1/2 ifK = Q{i) 

7/12 z/X = Q(y=3) 



l + \dK\ 



otherwise 



Proof For ^ just recall that t = VA/P by the definition of $ in (0) and 

t>tK. 

Thanks to the bounds on Re (z) and Im (z) given in the description of the 
fundamental domain J^k (in Definition 12. 5p we get 

• < |Rc(Q)| < P/2, < Im(-g) < 1/2, and so \Q\^ < /2 when 

• < Re(-Q) < P/2, -V3/6P < Im(-Q) < v^/3P and then \Q\^ < 
7/12P^ when K = Q(V^); 



• < Rc(-Q) < P/2, < Im(-g) < ^^p and then \Q\^ < 
In all cases we have 

IQI' < CkP^ < CK^2 

Recalling that PR — \Q\^ — A, we obtain 



P\ 



A 

K 



PR<[l + ^]l^. 



□ 



2.4 Julia's covariant 



From now on, let K be an imaginary quadratic field, let O be its ring of integers, 
and let Vq be the set of binary cubic forms in 0[x,y]. We want to define a 
canonical representative (or reduced form) in each orbit GL2(C') • P, P G Vq. 

Definition 2.11. We consider binary cubic forms in Vq, 

F{x, y) = ax^ + bx^y + cxy^ + dy^ , a,b,c,d O 
modulo the action o/GL2(0) given by 

M-F = {det{M))-^F{Ax+By, Cx+Dy), for each M = ^ e GL2{0). 

Remark. As we saw in CoroUarv 12.31 this is the restriction of the action used 
in Taniguchi's Theorem, when hpc — 1. 

Julia p^lj gives us a covariant for this action: 
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Definition 2.12. Let F G Vo be irreducible over K, factoring over C as 
F[x,y) = a[x — aiy){x — a2y){x — Oij,y), with a 7^ 0. We associate to F the 
positive definite binary Hermitian form 

HF{x,y) = tl\x - aiy\^ + tl\x - a2y\^ + tl\x - aaj/p, 

where 

t^ = IflPlo^j ~Q!feP, hiik pairwise distinct. 
The following three lemmas follow from a direct computation: 
Lemma 2.13. We have 

HF(x,y) ^P\x\^ + Qxy + Qxy + R\yf, 

where 

P^tl + t^ + tle R+, 
Q = -(ait? + ^2*2 + a3<i) e C, 
R^\ai\Hl + \a2\''tl + \a3\-'tleR+. 

Lemma 2.14. We have 

{t,t2t3)^ = \a\^\disc{F)\ (6) 
Lemma 2.15. Let A ^ - disc(i7F) ^ PR - |QP and D = disc(i^). Then 

A = 3|Z?|. (7) 
Proposition 2.16. The application which sends F to Hp is covariant, i.e. 

Hm-f — M ■ Hf, 

for allM e GL2{0). 

Thanks to this property we can translate our problem of defining a unique 
reduced F to the problem of finding a unique reduced covariant Hf plus some 
extra conditions as we will see in Section 13.21 

Definition 2.17 (Julia reduction). Let F — {a,b,c,d) E Vq be a binary cubic 
form with coefficients in O. We say that F is Julia-reduced (modulo GL2(C')J 
if its covariant Hp is reduced, in the sense of Definition \2.Si 

3 Reduction of binary cubic forms 
3.1 Bounds for binary cubic forms 

Let be a binary cubic form and let Hf be its covariant hermitian form. 
Starting from bounds on Hp coefficients it is possible to directly bound F 
coefficients and then to loop over all reduced binary cubic forms in time 0{X), 
but the coefficients involved in the complexity of this algorithm are quite big 
(see [13 for details) , so we chose another method, suggested by John Cremona, 
which can be found in [511^17]. 

Definition 3.1. For any k E O, we note Tk = Hi). 
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Definition 3.2. For any oq G O, we fix once and for all a system of represen- 
tatives Vaa for 0/3aoO. This is a finite set with 9 |ao|^ elements. 



Definition 3.3. Let Fk be as in definition \2.5\ 

We define Vk to he a fundamental region for C/O such that Fk C Vk- 

Proposition 3.4. Let F = {a,b,c,d) be a binary cubic form. There exists a 
unique fc G O such that sends F to an equivalent binary cubic form Fq = 
(ao, feoj cq, do) such that bo G Vao- We will call this Fq T-reduced. 

Moreover, if F is Julia-reduced, then we have also the following properties: 

\ao\ < 3-3/4^^3/2^1/4 

Icol < ^-^^ (8) 
either\do — Xi\ < ^^y= or \do — X2\ < ^ 



/\A\ ' " " - ^1-41' 

where ch = 3^/'^2-^/H]^^ , A = -27ag,B = ISaofooCo - 46g, C ^ b'^cl - Aoqc^, 
and we call Xi and X2 the roots of the quadratic polynomial Ax^ + Bx + C . 

Proof. As regards the first assertion, just remark tliat sends (a, b, c, d) to 
(ao, bo, Co, do) — (a, 6 + 3ak, 3ak^ + 2bk + c, ak^ + bk^ + ck + d). 

Now, assume tliat F is Julia reduced. 

Let us consider the seminvariants associated to Fq: 

Ph = b^ — SfloCo and Ur = 26q + 27agdo ^ 9aofeoCo 

(recall that Ph is the first coefficient of the Hessian of Fq, but it is not in general 
equal to Po, the first coefficient of the covariant associated to Fq). 

Tfc leaves Ph and Uh unchanged and, as shown in Womack's thesis [H] we 
have 

|ao| < 3-3/4t^3/2^i/8 (g) 



and 

so from the syzygy 
we obtain 



'^Ph = C/ff + 27disc(Fo)a^ 



Ph < chX^/\ (10) 

where ch — 'i^^'^2~^l^t'^ , and we easily obtain the bound for |co|. 

Finally, since disc(Fo) = A(f^ + Bdo + C and |disc(Fo)| < \fX we have 

\do-xr\\do-X2\ < VX/\A\, 



and this inequality implies that [do — ccij and [do ^ 2;2| can't be both bigger than 

□ 



Corollary 3.5. It is possible to list all the reduced binary cubic forms (a, b, c, d) 
(modulo GL2(0);, with Af{disc{F)) < X in time 0{X^+''), for all s > 0. 
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Proof. The number of r-reduced binary cubic forms (ao, 607 cq, do) which are 
equivalent to Julia-reduced ones (i.e. satisfying all properties enumerated in 
Proposition 13. 4p is 



Thus 



And 



^« E E E El 

or|do-an2|<SXi/''/|ao| 



^« E \<'o\'-—^-—^ = x- ^ 



|ao|«Xi/8 l"Ol n=l 

Since #{ao G O: |ao|^ = n} = O(n^) = 0(X'^) for all £ > 0, and J2n=i ^ is 
0(log(D)), we can conclude. □ 

Proposition 3.6. Let F be a Julia-reduced binary cubic form, let Fq be the cor- 
responding T -reduced form, and let Hp^ — (Pq^Qq, Rq) be the binary Hermitian 
form associated to Fq. Then F = T-k ' ^O; for a unique k € O. 

Proof. The action of sends Hp ~ (P, Q, R) to Hpg — {Pq, Qq, Rq) such that 
Pq = P and Qo = Q + kP. Dividing by P, we obtain zq = z ~ k, with 
z G Fk C Vk, but this uniquely determines fc, so we can conclude. □ 

Algorithm 3.7 (r-reduction). Let K be an imaginary quadratic number field of 
class number 1. This algorithm loops over all t -reduced binary cubic forms Fq ~ 
(flo, &0j Co, do) satisfying the conditions in Proposition \3.J\ with A/"(disc(Po)) < X , 
and associates the equivalent binary cubic form F — (a, b, c, d) = r_fe • Po; such 
that z G Fk (as explained in Proposition[ 



For each aQ,bQ,CQ,do in O satisfying the following properties: 

/ \ 3/2 

• 60 belongs to Vao 

• i^oi < — sra — , 

• Either |do - xi\ < X^/^y^] or \do - X2\ < X^l^l^\ 
Do the following operations: 

1. compute the first two coefficients Po,Qo of the covariant FLp^^ of the cubic 
form Po = (ao, 60, Co, do)- 

2. Compute k such that zo + G Vk {zq = —Qq/Pq). 

3. Compute P = (a, 6, c, d) — T_fc(ao, 60, co, do). 
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3.2 Automorphism matrices 

In this section we are going to study automorphism matrices for binary hermi- 
tian forms. 

Proposition 3.8. Let F = (a, 6,c, d), F Julia-reduced. Let H = Hp, and 
A = Pi?-|Q|'. 

Let M — (^'^ ^ ^ ^ G\j2{P) such that M ■ H = H . Then we have the 
following bounds on the coefficients of AI: 



and 

^ \B\<2^ ifB = Q 



B|<^ + 1 ^fB^Q ^^2) 



Proof. Let us write H{x, y) — P\x\'^ + Qxy + Qxy + R\y\^. We have 

PH{x,y)^\xP + yQ\^ + A\y\^, and (13) 

RH{x,y) = \Ry + Qx\^ + A\xf. (14) 

Thanks to formula ([T^ we can give upper bounds for and \D\. Let us 

write more exphcitly the relation M ■ H ^ H: 



I ApP + ACQ + ACQ + ICpfl ABP + ADQ + BCQ + CDR 
ABP + CBQ + ADQ + CDR \B\^P + BDQ + BDQ + \D\'^R 

H{A,C) 



H{B,D) 

By imposing this matrix to be equal to M we have 

l>ip + cgp + A|c|2 = p2 ^ \c\<JL^ 

V A 



P7? 

|BP + D0|2 + A|i?|2 = PP ^ |i5|2<^^ 

PR 

|CP + AQ|2 + A|^|2 = Pi? ^ < 

When C = the third equation becomes 

ABP + ADQ = Q, 

with \A\ = \D\ = 1 so it is easy to check that \AC\ P < 2Q < 2^fcKP and we 
obtain the formula. Finally, when C ^ 0, since \AD — BC\ = 1 we get 



\C\ 

and we easily conclude. □ 
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The bounds of the previous Proposition are completely explicit when hn = 1, 
since we know Ik and ck- 

Definition 3.9. Let M e PGL2(C'). We define 

S{M) = {H e ^/M; \M -H = H andH reduced}. 

that is, the set of reduced binary hermitian forms which are stabilized by the 
action of M . 

The following algorithm lists the finite set of automorphism matrices. It needs to 
be run only once for each of our 9 imaginary quadratic fields of class number 1. 

Algorithm 3.10. Computes the set Ai of all matrices M stabilizing some re- 
duced binary hermitian form, and for each M outputs also the corresponding set 
S{M). 



Set X = 0. 

For each triple {A,C,D) satisfying the bounds of Proposition 
lowing operations: 



do the fol- 



1. For each B e O such that \AD - BC\ = 1 (if C = 0, take only the set 

A B 



\B\ < 2y/c^), let M = ^ ^ j and do the following. 
2. Consider the following 4x4 matrix, with coefficients in O: 



W{M) 



( iW-i) 

AB 
AB 
V \B\' 



_AC 

{AD- 1) 
CB 
BD 



Aq_ 

BC 
{AD- 1) 
BD 



CD 
CD 



3. Compute the rank r of W{M) (over the field K). 

4. If r = 1 or r = 4, skip to the following quadruple {A, B, C, D). 



5. If r = output M = 

H reduced}. Myj{M} 



A B 
C D 



and S{M) ^ {H e 



6. If r = 2 or r = 3, set M = 

{P.Q.R) e 

S{M)t^%, output M = 



A B 
C D 

X \ W ■ {P,Q,Q,RY = and {P,Q,R) reduced}. If 
A B 



C D 



compute the set S{M) = {H 
= and {P, Q, R) reducei 
and S{M). M=MU {M}. 



Output M. 
Remarks. 

• We could also loop only on A, D and replace step (1) by : 

1. Solve \AD - BC\ = 1 for B, C € O. This time BC belongs to an 
explicit finite set, and we enumerate divisors. 
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• It is possible to write (once for all) explicit conditions to associate to any 
binary hermitian form H its set of automorphism matrices, just looking 
at the sets S{M) computed in algorithm 13. 101 

• For an example of application of the above algorithm, Appendix Elcontains 
the list of all automorphism matrices for K = and the corresponding 
conditions on binary hermitian forms. 



Remark. Running the algorithm on all the 9 possible number fields we noticed 
a property holding for all K ^ Q(\/"-3) : 

• For each matrix M found at step 6 (that is, they are not trivial automor- 
phisms) VF(M) has rank 2 and S{M) is a subset of the boundary of the 
fundamental domain. 

In the case K = Q(V— 3) we have explicit counter-examples. 



The proof of Algorithm 13. 101 is given by the following proposition. 

Proposition 3.11. Let M = (^ |^) g PGL2(0) belong to the stabilizer of Hp, 
where Hp is the Hessian of some reduced cubic form F. If r is the rank of the 
matrix W constructed in the above algorithm, then 

• r — if and only if B = C = and A — D are units. Then M is an 
automorphism for all Hermitian quadratic forms in T . 

• r = \ is impossible 

• r — 2 or r = 3 then M is an automorphism for some linear subspace of 
^ , defined by explicit equations in the variables P,Q,Q, R. 

• r = A is impossible. 

Proof. The condition {cd)^ Aut(iJ) translates to the linear system W{M) ■ 
X = Q, with X = {P,Q,Q,RY. 

• If r = 4, the only solution of the system is (0,0,0,0), but this is not 
allowed since P,R>0. 

• Assume that r < 1 : the matrix ( ^ £ ) has rank 2 so the two 2 by 2 
matrices on the lower- left and upper-right corners of W{M) have rank 2 
unless B = C = 0. In this case W{M) is diagonal 

/ |Ap-l _ \ 
AD-1 

AD -I 

Since i? = C = 0, and AD — BC is a unit, we must have \A\ = \D\ — 1, 
so this matrix has either rank 2 or (when AD ~ AD ~ 1). 

□ 
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4 The algorithm 



Algorithm 4.1. Given a hound X — , output the list of reduced binary cubic 
forms modulo GL2(0), such that Af {di8c{F)) < X. 

Use sub-Algorithm im to loop over quadruples F — {a,b,c,d) € satisfying all 
the properties in Section \JA[ Do the following operations: 

1. Approximate the complex roots of F, {ai,a2,a3) to a sufficient accuracy. 
Then approximate H = Hp = {P,Q,R) the associated Hermitian form. 

2. Check if H is in the fundamental domain modulo PGL2(C') (i.e. it is re- 
duced), (see Definition 12. Sp . In particular, \f Hp is "near" to the boundary of 
the fundamental domain use Algorithm 15.2! (see below) to check exactly the 
boundary condition. If not skip to the following F. 

3. Check whether F is irreducible in K[x,y]. If not skip to the following F. 

4. Apply Dedekind criterion to check whether F describes a maximal ring. If 
not skip to the following F. 

5. Apply sub-algorithm l3.10l to compute A^, the set of all automorphism matrices 
for H. 

6. Compute the set {M ■ F \ M G Ai} and check if F is the minimal element 
of this set (for some order, for instance the lexicographic one). If not skip 
to the following F. 

7. print F. 
Remarks. 

• For the precision needed in step (1) refer to Appendix C of [15] . 

• In step (5), we compute a list of automorplrs for F to decide wlietlier 
F is minimal among the reduced forms in its orbit with respect to the 
lexicographic order (in this case F should be kept, otherwise not). Another 
way to deal with this problem would be to store all those F and then check 
GL2(0)-equivalence once we have all the forms with a fixed discriminant 
D. The problem is that our algorithm does not output forms ordered by 
discriminant, so we could apply this test only at the end, and this would 
increase dramatically the space complexity. (Remember that we output 
the "good" binary cubic forms as we find them, so we do not keep in 
memory the list of representatives of cubic extensions) . 

5 Implementation problems 

5.1 Checking rigorously the boundary conditions 

As the computation of P,Q,R involves floating point approximations of the 
complex roots of a polynomial in ©[A], it will not give, of course, exact results. 
Those floating point computations will in general be sufficient to test whether 
the Hermitian form is strictly inside or outside the fundamental domain. But if 
it is very near the boundary (or worse on the boundary), this approach fails. 



14 



To get rid of this problem we use the fohowing formulas: 



\b\' 



P = -j^ + 3(|ai|^ + |a2|^ + |a3r) (16) 

be 

Q = -j— 75- + 3(Q;ia2a3 + aia2a3 + aia2a3) (17) 



c 



2 



R = -^ + 3(|ai|^|a2r + |air|a3r + |a2rl«3r) (18) 

Now we consider ai, 02, a3, cET, 02, 03 as algebraic numbers, and we let S be 
the set of the six permutations fixing the a^, and acting as on the al. The 
polynomial 

gp =W{^^ cr(aiaT + 0:202 + aacEs)) 
cres 

vanishes at |aip + |Q;2p + |a3p, and its coefficients are symmetric in (ai, a2, 0:3) 
and (oT, 02 , 03) independently. They can thus be expressed in terms of (6/a, c/a, d/ a) 
and (b/a, c/a, d/a). The polynomial fp{X) — gp — 3|l!p) vanishes at P and 
belongs to K[X]. 

In the same way we can compute polynomials in iir[^] vanishing at Q, 
R, Re((5) or Im(Q). Such polynomials are easily computed using a computer 
algebra system like Maple (and it is sufficient to compute them once for all). 

We want to verify rigorously boundary conditions, for instance P = R: if fp 
and fp have no common factor in then P ^ R. But this is not enough: 

we also want to check whether P < R or P > R, i.e. if the point we are testing 
is "inside" or "outside" the fundamental domain. 

The following theorem of Mahler [12J provides the accuracy we need for our 
floating point computations: 

Theorem 5.1 (Mahler). Let f — qqx"'' + aix"^^^ + ■ ■ — ao{x — ai) ■ ■ ■ {x~ 

am) be a separable polynomial of degree m > 2, and let 

^(/) = , min \ai-aj\ 

be the minimal distance between two distinct roots of f . Then 
A(/) > V3m-(™+2)/2|disc(/)|i/2M(/)-(™-i), 

where disc(/) is the discriminant of f , and M{f) = |ao| 11™=! ™^-^(-'^' I'-'^'il)- 
This translates to the following algorithm: 

Algorithm 5.2 (Checking an algebraic identity). Let a and /3 e M 6e two 
algebraic numbers, and let A and B £ K[X] \ that vanish at a, and /3 respec- 
tively. Assume we can compute floating point approximations a and (3 such that 



a\ < e, 



(3-(3 



< e, for any fixed £ > 0. 



We want to decide whether a<l3,a>l3ora — p. 

1. Let C = AB snd f = C/gcd(C, C"). 

2. If the degree of / is 1, then answer a = (5. 
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3. Compute a good approximation A of 

A(/) = V3m-("+2)/2|disc(/)|i/2M(/)-(™-i), 
where disc(/) and M{f) are defined in Theorem 15.1! such that A < A(/). 

4. Compute a and /3 at precision e ~ A/4, i.e. a and j3 such that 

/3-/3 



la — a| < e, 



< £. 



5. If jd — /3| < 2e, answer a — (3. 

6. If ci < (3, answer a < (3. 

7. If d > /?, answer a > /?. 

Proof. The polynomial / is non constant and has a and /3 among its roots. If 
its degree is 1, then a — (3. Otherwise, assume first that |d — /3| < 2e. Then 



/3|< 



P~f3 



a — P 



< 4e < A(/). 



Hence a = (3 hy Mahler's theorem in this case, proving 
We now assume that jd — /3| > 2e; since 



Q;-/3 = d- /3+(a-d)~(/3-/3) 



and 



a — (3 and a — (3 have the same sign. 



< 2e, 



□ 



Proposition 5.3. The smallest e that we can obtain in step (4) of the above 
algorithm (i.e. the maximal precision needed) is 3> X^^, for some positive 
constant (3. 

Remark. That means that for our computation we will need at most J7(logX) 
significant digits. 

Proof. Algorithm 14.11 loops over reduced integral cubic forms F — (a, 6, c, d) G 
Vo with discriminant disc(F) satisfying A/'(disc(i^)) < X. In particular, Propo- 
sition [S3] implies that \a\ < X'^l^. 

For each such form, we may compute various separable polynomials / with 
coefficients in a^^Ox, for some bounded integer u. Then disc(/) is non zero, 
in o^'^'^Ok- Its norm is a non-zero rational integer divided by \a\ hence 
> X-'". Thus disc(/) > X-"/^ 

Landau's theorem (see [3, Proof of Theorem 13.1] for example) tells us that 

Mif) < II/II2 

and the coefficients of / are monomials in ei, 62, 63, /i, /2, /s (see Appendix D 
of [TS]). Each one of these is bounded by c • AT", for an appropriate constant c 
and exponent a. 
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We have 

A(/)» Af(/)-(™-i). 

So we obtain 

ll/l|2«^^ 

but then we can conclude that A(/) > X^^ . □ 

5.2 An idea to count only half of the extensions 

Let K be an imaginary quadratic number field, with class number hx = 1 and 
discriminant cLk ^ —3, —4. It is easy to remark that if H — (P, Q, R) is in the 
fundamental domain, then also H' = {P, —Q, R) is. And, in general, these two 
Hermitian forms are not equivalent modulo PGL2(C'). 

In particular, ii F = {a,b,c,d) has Hp = H, then F' = {a, —b,c, —d) gives 
Hf' = H'. 

So we can loop only on half of the c satisfying the given bounds, then con- 
struct both the forms F = (a, 6, c, d) and F' = (a, — 6, c, —d) and check if they 
are equivalent (comparing F' with the list of automorphic functions to F) . If 
not we verify also the list of automorphic functions to F' to see if one of them 
will be found in our loops, and if both answers are no, we add this second form 
F' to our output list. 

6 Results 

In this section we present results obtained for the case K = Q(i) via an im- 
plementation of our algorithm in Pari/GP [16J, running on an Intel Xeon 5160 
dual core, 3.0 GHz. 

Let X be the bound on M {X){L / K)) and N{X) the number of isomorphism 
classes of cubic extensions of Q(i) up to that bound. 

In the following table, we will compare the time needed to list all N{X) 
representatives of cubic extensions of Q(i) with two algorithms : ray class fields 
one and ours, t denotes the running time of our algorithm, t' the running time 
of ray class field one (see Section 9.2.3 of [3]). 



X 


N{X) 


t 


t' 


10* 


276 


5 s 


16 s 


4- 10* 


1339 


19 s 


Imn 18 s 


9 • W 


3305 


56 s 


3mn 45 s 


10*^ 


42692 


24 mn 1 s 


2h 52mn 9 s 


4- 10« 


181944 


2 h 49 mn 


34h 24 mn 8 s 


9 • 10« 


421559 


9 h 37 mn 


> 134 h 


10« 


4990974 


359 h 25 mn 


> 2720 h 



Remarks. 

• This computations allowed us to check the correctness of our results. In 
fact, we compared N{X) up to X — 9-10^ with the results of the ray class 
field algorithm and all results matched. 
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• The last line of the table would have involved too long computations with 
ray class field algorithm, so we skipped it, and we give only a prevision on 
the running time needed. 



6.1 Roberts' conjecture and asymptotic predictions for 

N{X) 

Frank Thorne compared our numerical results with heuristic asymptotic de- 
velopments derived from the Datskovsky- Wright method [8J, in the spirit of 
Roberts's conjecture (see [E]). Starting from Roberts' conjecture, Taniguchi 
and Thorne worked out the formula in the particular case when k is an imagi- 
nary quadratic number field : 



1 ReSs=iCK{s) ^ , 1 ,-1/2^^^ , , , V3r(l/3)6 Ck(1/3) 



The following table compares our values for N{X) with Thome's asymptotic 
data. The results are strikingly similar. 



X 


N{X) (Morra) 


N{X) (Taniguchi-Thorne) 


10* 


276 


270.2 


10*^ 


42692 


42655.6 


9 • 10" 


421559 


421260 


10« 


4990974 


4990962 
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A Automorphism matrices for Q(i) 



M (modulo multiplication by £ e {+1, — — «}) 



conditions for S{M) 



1 
1 



for all P,Q,R 



-1 

1 



iiP = R and Re(g) = 



1 

1 



if P = i? and Im(Q) = 



1 

-i 



iiP = R and Re(Q) = Im(Q) 



1 

1 



].iP = R and Re(Q) = - Im(Q) 



-1 

1 1 



1 1 
-1 



if P = i? and Re(Q) = P/2 



-1 

1 -1 



-1 1 

-1 



iiP = R and Re(Q) = -P/2 



1 

1 / 



-i 1 
1 



if P = P and Im((7) = P/2 



-1 
1 1 



if P = P, Re(g) = P/2 and Im(Q) = 



i W -1 -l-i 
-i 1 j'l 1 1 
-i \ / i 

1 1 y ' V 1 + 1 

1 -i + 1 
-1 



1 1 

-1 



1 -i 
-1 



-1 -\-% 

-i 1 
1 1 
— z — z — 1 
1 
i-1 -1 
1 

-1 i-1 



if P = P,Re(Q) = P/2 
and Im(g) = P/2 



-1 

-1 1 



if P = P, Re(Q) = -P/2 and Im(Q) = 



-i 
-i 1 
i 
-1 1 

1 

1 ' -i 
l-i -1 



-1 l-i 
-1 1 

-i 

1 -i-1 
-1 - i 

-1 

1 -1 
l-i -1 



-1 l-i 
-i 1 
1 -i 
1 -i-1 
1 
l-i -1 

1 -1 
—i i — 1 



if P = P,Re(g) = -P/2 
and Im(g) = P/2 



-1 

-i 1 



if P = P, Re(g) = and lm{Q) = P/2 





-1 



1 
-i 



1 
i 



if g = 



1 1 

-1 



if Re(g) = P/2 and Im(Q) = 



1 i 

-i 



1 1 + i 
-1 



if Re(g) = P/2 and Im(g) = P/2 



if Re(Q) = -P/2 and Im(Q) = 



1 i 
i 



-1 
—i 



1 -1 




+ 1 
1 



if Re(Q) = -P/2 and lm(0) = P/2 



1 i 
-1 



if Re(g) = and lm{Q) = P/2 
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